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ABSTRACT 

We present results on gas flows in the halo of a Milky Way-like galaxy at z = 0.413 based on high- 
resolution spectroscopy of a background galaxy. This is the first study of circumgalactic gas at high 
spectral resolution towards an extended background source (i.e., a galaxy rather than a quasar). Using 
longslit spectroscopy of the foreground galaxy, we observe spatially extended Ha emission with circular 
rotation velocity Vdrc ~ 270 km s“^. Using echelle spectroscopy of the background galaxy, we detect 
Mg II and Fe ii absorption lines at impact parameter p = 27 kpc that are blueshifted from systemic 
in the sense of the foreground galaxy’s rotation. The strongest absorber (EW 2796 = 0.90 A) has an 
estimated column density {N^ > 10^® cm“^) and line-of-sight velocity dispersion (u = 17 km s“^) 
that are consistent with the observed properties of extended H i disks in the local universe. Our 
analysis of the rotation curve also suggests that this r « 30 kpc gaseous disk is warped with respect 
to the stellar disk. In addition, we detect two weak Mg ii absorbers in the halo with small velocity 
dispersions (cr < 10 km s“^). While the exact geometry is unclear, one component is consistent with 
an extraplanar gas cloud near the disk-halo interface that is co-rotating with the disk, and the other 
is consistent with a tidal feature similar to the Magellanic Stream. We can place lower limits on the 
cloud sizes {I > 0.4 kpc) for these absorbers given the extended nature of the background source. 

We discuss the implications of these results for models of the geometry and kinematics of gas in the 
circumgalactic medium. 

Subject headings: galaxies: evolution - galaxies: halos - galaxies: ISM 


1. INTRODUGTION 

The interplay between inflows and outflows of gas 
around galaxies has wide-ranging implications for galaxy 
evolution. Of particular interest are the roles of gas in- 
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Cool gas clouds (T ~ 10^ K) in galaxy halos are of par¬ 
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individual Mg ii absorbers with individu al galaxies (e.g., 


Chen et al. 2010b Nielsen et al. 20131, but it is clear 
that the average equivalent width and covering fraction 
are larger at smaller impact parameters. In particu¬ 
lar, Mg II halos extend to p « 100 kpc for L* galax¬ 
ies with a covering fraction of » 50% for absorbers with 


V1V(2796) > 0.5 A (Chen et al. 


fraction at p < 25 kpc approac 

eTalllM^ . 


2010a), and th e covering 


les unity (e.g., Kacprzak 


Beyond these trends with impact parameter, there is 
strong evidence that the equivalent width and cover¬ 
ing fraction of Mg ii absorbers depends on the orien¬ 
tation of the background sightline with respect to the 
disk of the foreground galaxy. For example, when stack¬ 
ing low-resolution spectra of Zr,hnt > 1 ga l axies behind 


0.5 < Zspec < 0.9 galaxies, Bordoloi et al. (2011) found 


that edge-on disk galaxies at small impact parameter 
(p < 40 kpc) exhibit three times stronger absorption 
along the minor axis than along the major axis. This is 
consistent with the signature of bipolar outflows driven 
perpendicular to the plane of the disk, as observed in 
local starburst galaxies (e.g.. Bland fc Tully||1988 Shop- 


bell & Bland-Hawthorn 199S'Uhyama et al.| 20l)2|). It 
is also consistent with the incidence of outhowmg gas 
traced by Mg ii absorption in down-the-barrel sightlines 
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motivated studies of the spatial extent, mass, and kine¬ 
matics of outflowing gas for individual galaxies based on 


sightlines towards background quasars (e.g., Tripp et al. 
20111 IBouche et al.||2012[ [Gauthier fc Chen]|2012[ |Lund- 
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In addition to a preference for absorption along the mi¬ 
nor axis, there is evidence for Mg ii absorpt ion associated 
with t he p lane of the disk. In pa rticular, [Bouche et al.| 
(2012) and Kacprzak et al. (2012|) argued that the distri- 
bution of azimuthal angles for individual Mg ii absorbers 
is bimodal, with gas being found preferentially within 
« 20° of either the minor or major axis. While it is chal¬ 
lenging to characterize the detailed three-dimensional ge¬ 
ometry of cool halo gas based on individual sightlines, 
there are clear examples for which Mg ii absorbers are 

co-rotating with the stellar disk (e.g., jSteidel et al.|2002| 
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al.| |2013p . For example, |Steidel et al.| (|2002p compared 

the kinematics of Mg ii absorbers at p = 10-50 kpc to 
emission-line rotation curves at r < 10 kpc and found 
that the absorption-line velocities were often offset from 
the systemic redshift in the sense of the galaxy’s rotation. 
In detail, the full velocity ranges of Mg ii absorption 
(Au « 100-200 km s“^) were broader than expectations 
for a single thick disk, but the overall kinematics were 
clearly dominated by rotation. 

While absorption-line studies of halo gas have typi¬ 
cally used sightlines to background quasars, the advent of 
large redshift surveys and sensitive spectroscopy with 8- 
10 m class telescopes has enabled the use of background 
galaxies as li ght beacons. This has been done in both 
stacked (e.g., Steidel et al. 2010 Bordoloi et al. 2011) 


and individual s 

oectra (e.g., 

Adelberger et al. 112005 

Ru-| 

bin et al. 

2010a 

Steidel et ai.||2010) at low spectral reso- 


light profiles of background galaxies provide an opportu¬ 
nity to study absor ption over a larger s olid angle through 
a given halo (e.g., Rubin et al.||2010a), and the number 


density of galaxies on the sky otters significantly more 
sightlines to study halo gas in a bsorption at faint mag¬ 
nitudes (e.g., Barger et al.||2008 ). 


However, it has not yet been possible to obtain spectro¬ 
scopic data for individual background galaxies that are 
of comparable quality (e.g., S/N, spectral resolution) to 
what is available for background quasars. In particular, 
high-resolution spectroscopy is required to measure the 
kinematics, column density, and covering fractions for 
individual absorption components, and these measure¬ 
ments provide important constraints for models of cool 
gas clouds in the circumgalactic medium. 

In this paper, we use a bright galaxy from the 


Diamond-Stanic et al. (2012) sample as a background 
light source [Zbg = U.712) to probe halo gas associated 
with a disk galaxy (zfg = 0.413) at impact parameter 
p = 27 kpc. The brightness of the background galaxy 
{g = 19.6) offers the first opportunity to study halo gas 
at high spectral resolution towards an extended back¬ 
ground source. We describe the data used in this study 
in Section 2. We present results on the properties of the 
foreground disk galaxy and its circumgalactic gas in Sec¬ 
tion 3 and 4. We discuss the implications of these results 
in Section 4, and we summarize our main conclusions in 
Section 5. Throughout the paper, we adopt a cosmology 
with h = 0.7, Dm = 0.3, and Da = 0.7. 

2. DATA 

2.1. HST/WFC3 

Observations with the Wide Field Camera 3 on the 
Hubble Space Telescope were obtained on 20 10 Dec 15 
as part of program 12272. As described by |Diamon(T| 
Stanic et al. (2012), we obtained 4 x 10 min exposures 
in a single orbit using the F814W filter on the UVIS 
channel (pixel size 0.04"), and we combined the dithered 
images to produce a science mosaic with 0.02" pixels. An 
11" X 11" portion of the image including the foreground 
disk galaxy and the compact background galaxy is shown 
in Figure [2 

2.2. Keck/LRIS 

Observations with the Low R esolution Imaging Spec¬ 
trometer (LRIS,|Oke et al.|1995) on the Keck I telescope 
were obtained ori 2011 iVlar 7. Using a 1" longslit and the 
D560 dichroic, we obtained a 30-min exposure with the 
400/3400 grism on the blue side (i? « 600) and 2 x 15 min 
exposures with the 400/8500 grating on the red side {R « 
1000). The slit position angle was set to -41.9 deg to in¬ 
clude both the background (SDSS J090523.59-1-575912.4) 
and foreground (SDSS J090524.08-1-575915.9) galaxies 
(see the solid red lines in Figure 1). The data were pro¬ 
cessed with the XIDL LowRedu> ^ pipeline. 

2.3. Keck/HIRES 

^ http://www.ucolick.org/~xavier/LowRedux/ 
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Fig. 1. — HST WFC3/UVIS F814W image of the background and foreground galaxies described in this paper. The positions and 
orientations of the Keck/LRIS (solid red lines) and Keck/HIRES (dashed blue rectangle) slits are shown. The orientation of the foreground 
disk is shown by an ellipse, and the direction of rotation is indicated by the arrow on the near side of the disk. The dotted portion of the 
ellipse corresponds to the far side of the disk, which is oriented into the page. 


Observations of the background galaxy w ith the High 


Resol ution Echelle Spectrometer (HIRES, Vogt et al. 


19941 on the Keck I telescope were obtained on 2U(J9 Dec 


TE We obtained 2 x 1 hr exposures using the HIRESb 
conhguration with a 1.148" x 7.0" slit (see the dashed 
blue rectangle in Figure [R, which provides a spectral 
resolution of i? « 37, OOC) and a velocity resolution of 
8 km s“^. The data were processed with the XIDL 
HIRedujfj pipeline. The spectrum of the background 
galax y is lairly smooth near the f oreground absorption 
lines (Diamond-Stanic et al. 2012), and we fit the con- 
tinuum tor each order using a low-order polynomial, ex- 


http://www.ucolick.org/'xavier/HIRedux/ 


eluding regions of line absorption and emission. The 
continuum-normalized data from the two exposures and 
from overlapping regions between the orders were com¬ 
bined using an inverse-variance weighted mean. 


2.4. Spitzer/IRAC 
Observations w ith the Infrared Array Camera (IRAC, 


Fazio et al. 2004) on the Spitzer Space Telescope were 
obtained on 2009 Nov 27 as part of Program 60145. We 
obtained 5 x 30 sec dithered exposures in the 3.6 pm 
and 4.5 pm bands. We performed photometry on the 
post-basic calibrat ed using the MOPEX point-s ource ex¬ 
traction software (Makovoz & Marleau 2005). Given 
the close spatial separation of the background and fore- 
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Fig. 2. — Left: This color-magnitude diagram shows the location of the foreground galaxy with respect to a comparison sample of 
~ 50, 000 galaxies at comparable redshifts from the PRIMUS survey. The contours encompass 50%, 75%, and 90% of the galaxies from the 
comparison sample. The U — B color of the foreground galaxy place it near the luminous end of the green valley. Right: The star-formation 
rate vs stellar mass plane for the same set of galaxies, which illustrates the location of the foreground galaxy near the massive end of the 
star-forming sequence. 


ground galaxies (5.2") and the spatial resolution of IRAC 
(FWHM = 1.7"), we used point-source photometry to 
properly deblend the two galaxies. 

3. PROPERTIES OF THE FOREGROUND DISK 
GALAXY 

3.1. Redshift, Rest-frame Color, Stellar Mass 

We measure a systemic redshift of z = 0.4134 ± 0.0001 
for the foreground disk galaxy based on [O ii] A3727, 
[O III] A5007, and Ha A6563 emission lines in the 


Keck/LRIS spe ctrum. We then use K-correct (Blanton 
& Roweis 2007) to estimate rest-frame color {Jj — B = 
1.01) and absolute magnitude {Mg = —20.91) based on 
photometry from the Sloan Digital Sky Survey (SDSS 
Abazajian et al. 2009) and from Spitzer/IRAC. We es¬ 
timate the stellar mass (Al* « 2.0^°"^® Al©) and star 
formation rate (SFR « 2.4 Al© yr“^) from the same 
photometry using iSEDfit, the spectral energy d i stribu - 
tion modeling code described by|Moustakas et al.|(|2013|). 
Thes e values are derived fo llowing the hducial method 
from Moustakas et al. (2013), which uses Flexible Stellar 
Population Syn th esis models ( Gon roy et al. 2009 Conroy 
& Gunn 2010), a Chabrier (2003) initial mass function 
from 0.1 to 1110 A l^. ancUtne time-dep endent dust at¬ 
tenuation curve of|Charlot & Fall] (|2000D. 

We show the location ot the foreground disk galaxy on 
a color-magnitude diagram and a star formation rate- 
stellar mass diagram in Figure As a comparison sam¬ 
ple, we include « 50,000 galaxies at 0.3 < z <0.6 from 
the Prism Multi-object Surve^FI (PRIMUS, ICo il et al. 
2011at ICool et al.l 20131 llVloustaKas et al.l 20131). This 


with a star-forming galaxy viewed at moderate inclina¬ 
tion. The stellar mass and luminosity of the foreground 
galaxy place it near M* for the stellar mass function 


(Moustakas et ajj 

2013) and L 

for the luminosity func- 

tion at z ^ U.4 ( 


Gilmer et al. 

|2006). In particular, the 
h i l' n QAvy r* 


by the H-band luminos ity function fo r (blue , red, all) 


galaxies at z = 0.5 from Willmer et al. (2006). 


3.2. Rotation curve and dynamical mass 

We show the two-dimensional Keck/LRIS spectrum in 
the region surrounding the Ha emission line for the fore¬ 
ground galaxy in Figure ^ The Ha emission is spatially 
extended and exhibits the signature of rotation with a 
line-of-sight amplitude vios = 178 ± 10 km s“^. We can 
convert this into a circular rotation velocity (vrot) by 
accounting for the disk inclination (i) and the azimuthal 
angle of t he slit relative to th e m ajor axis of the disk (a). 
Following Chen et al. ( 20051) and Chen et al. (2014), this 
conversion can be expressed in terms ot in terms or i and 
a for any location in the disk plane. 


'^rot vT -|- sin^ a tan^ i 


vios cos a sm i 


( 1 ) 


We fit an exponential pro file to the HST data with 
GALFIT (Peng et al. 2002) to estimate the axis ratio 
{b/a = 0.39) and position angle of the disk (PA = 15.8°). 
This implies a disk inclination i = 67° and an azimuthal 
angle a = 33° of the background sightline relative to the 
major axis of the foreground disk. Given the system¬ 
atic uncertainties associated with e stimating projection 


hgure illustrates that the foreground disk galax y resides param eters for spiral galaxies (e.g., Barnes & Sellwood 


et al. 2011) and the massive end of the star-forming se- 
quence (e.g., Noeske et al.||2007 ). The dust attenuation 
we estimate from iSEDht is Ay « 0.5, which is consistent 

http://primus.ucsd.edu/ 


on the lurn inous end of the green valley (e.g., Mendez 2003), we adopt ±5° uncertainties ori the inclination and 


azimuthal angles. 

Before using equation [J to estimate a circular veloc¬ 
ity, we need to account for the fact that each spatial 
bin of the LRIS spectrum includes emission from re¬ 
gions that span a large range in azimuthal angle. For 
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Fig. 3.— The two-dimensional spectrum from Keck/LRIS near AqI^^ = 9276 A, showing the Ha emission from the foreground galaxy 
(top) and the bright continuum of the background galaxy (bottom). The Ha emission from the foreground galaxy is spatially extended 
and exhibits the signature of rotation. The ionized gas in the disk of foreground galaxy is blueshifted by u ^ —200 km s“^ on the side 
closest to the background galaxy. 


example, Figure H shows that the brightest regions of 
the galaxy that fml within the slit are close to the ma¬ 
jor axis (i.e., they have small a angles). To account for 
this, we determine all the locations in the galaxy that 
would fall in each 1.0” wide x 0.2” long spatial bin along 
the LRIS slit. We calculate the relevant a angle for all 
500 pixels in the bin and then determine the appropri¬ 
ate Vrotivios ratio by computing a flux-weighted aver¬ 
age, assuming that the Ha morphology is similar to the 
observed HST morphology at Arest ~ 5800 A. For the 
spatial bins that correspond to the flat part of the ob¬ 
served rotation curve, we find that this ratio converges to 
Vrotjvios = 1.51 ± 0.05, which corresponds to a circular 
rotation velocity v^ot = 269 ± 18 km s“^. 

The GALFIT prohle fitting also yields an estimate of 
the half-light radius, = 5.7 ± 0.3 kpc. We estimate 
the uncertainty in this value from variations between fits 
that include Sersic profiles (n is a free parameter) and ex¬ 
ponential profiles (n = 1), and fits that include a second 
component to model emission from the central kilopar- 
sec of the galaxy. For the two-component fits, the fainter 
central component contributes 7 ± 2% of the total light 
and has re = 0.4 ±0.1 kpc. The best-fit Sersic index for 
this central component is n = 0.6 ± 0.1, which indicates 
that it does not exhibit the characteristics of a classical 
bulge (i.e., it is not well fit by an n = 4 de Vaucouleurs 
profile). 

We note that a circular velocity Vc = 269 km s“^ at 
Te = 5.7 kpc corresponds to a dynamical mass Mj_yn(r < 


Te) « TeV^jG = 9.6 X 10^° Al©, which is 0.5 dex larger 
than our estimate of the stellar mass above. This offset 
is somewhat larger than the 0.3 dex systematic uncer¬ 
tainty in the stellar mass estimate, so this suggests that 
the dynamical mass has comparable contributions from 
baryons and dark matter within the effective radius. 

4. PROPERTIES OF GAS IN THE 
GIRCUMGALAGTIG MEDIUM 

As discussed in the Introduction, the serendipitous su¬ 
perposition of the background and foreground galaxies 
on the plane of the sky presents an opportunity to study 
gas in the circumgalactic medium (GGM) of the fore¬ 
ground galaxy at impact parameter p = 27 kpc. In 
particular, we are able to measure the column density 
and kinematics of gas in absorption along the line of 
sight to the background galaxy at high spectral resolution 
(FWHM « 8 km s“^). In this section, we present these 
measurements, and we describe the kinematic modeling 
that is necessary to interpret the observed line-of-sight 
velocities in terms of a rotating disk, a warped disk, or 
an outflowing wind. 

4.1. Observed column densities and kinematics 

We show the continuum-normalized Keck/HIRES data 
for the Mg ii AA2796,2803 and Fe ii AA2586,2600 absorp¬ 
tion lines in Figure We also use spectral coverage of 
three other Fe ii lines (2344, 2374, 2382) in our analysis. 
We detect an optically thick absorption component cen- 
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Fig. 4.— Continuum-normalized spectra from Keck/HIRES showing absorption-line profiles for four transitions associated with the 
foreground galaxy (Mg II 2796, 2803; Fe II 2586, 2600). The solid red, green, and blue curves in each panel show the best-fit Voigt profiles 
for three components at u = —195 km s“^, v = —256 km s“^, and v = —363 km s"'^ relative to the systemic redshift of the foreground 
galaxy (z = 0.4134). The kinematics of the strongest component at u = —195 km s“^ suggest that it is consistent with an extension of a 
fiat rotation curve to r Ri 30 kpc. 


tered near v « —200 km s“^, along with two weaker com¬ 
ponents near v « —260 km and v « —360 km s“^. 
We use VPFITp^to perform a joint fit to all three compo¬ 
nents for both ions to measure precise column densities 
and kinematics (see Figure]^ and Table 0. In particular, 
we require that the redshifts (i.e., centroid velocities) for 
the three velocity component be the same for both Mg ii 
and Fe ii but we allow the b values and column densities 
to be fit independently for each ion. 

Using reasonable assumptions about gas-phase abun¬ 
dances and ionization, we can estimate the hydrogen col¬ 
umn density and [Mg/Fe] for each velocity component. 
For our fiducial estimates, we assume no dust depletion 

http://www.ast.cam.ac.uk/ rfc/vpfit.html 


and solar metallicity, and we assume that Mg ii and Fe ii 
are the dominant tracers of Mg and Fe, respectively. Re¬ 
laxing any of these assumptions would increase our es¬ 
timates of Nh (i.e., dust depletion or lower metallicity 
would reduce the gas-phase abundance of Mg and Fe rel¬ 
ative to hydrogen, and other ionization phases could con¬ 
tribute additional Mg and Fe). Based on the Fe ii column 
densities, these fiducial assumptions imply hydrogen col¬ 
umn densities Nh > 1-2 x 10^® cm“^ for the strongest 
component and Njj > 3-4x 10^^ cm“^ for the two weaker 
components. We note that the estimated dust depletion 
for Fe in the warm ionized medium is —0.8 dex (e.g., 
Jenkins 2009 Draine 2011[) and Lyman Limit Systems 


are typically associ ated with sub-solar metallicity (e.g.. 


Lehner et al. 2013), so these Nh values could be larger 
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by an order of magnitude. 

Based on the Mg ii and Fe ii column densities, all 
three velocity components have [Mg/Fe] « —0.2. The 
Mg II lines for the strongest component are saturated, 
which implies that there could be contributions from un¬ 
resolved velocity subcomponents. That said, the Fe ii 
lines span a much wider range in optical depth, and the 
b values we measure for the Mg ii and Fe ii components 
are nearly identical. This suggests that our column den¬ 
sity measurements for both ions are robust, even for this 
strongest component. We note that correcting for dust 
depletion would further dec rease [Mg/Fe] because Mg i s 
less depleted than Fe (e.g., Jenkins 2009 Draine 2011). 
The small [Mg/Fe] ratio for this absorbing gas suggests 
that it has been significantly enriched in Fe by previ¬ 
ous g enerations of type la supernovae (e.g., Rigby et al. 

loo^. 


4.2. Kinematics of gas in the disk plane 

If the kinematics of these absorbers were dominated 
by disk rotation in the plane defined by the stellar emis¬ 
sion, we could convert the observed line-of-sight velocities 
(vios) into deprojected rotation velocities (vrot)- This de¬ 
pends on the sine of the disk inclination angle (z) and the 
cosine of the azimuthal angle in the plane of the disk {(f). 


'^rot 

'^los 


1 


Sin i cos ( 


( 2 ) 


Following Steidel et al. (2002), the angle (j) can be ex¬ 
pressed in terms ot the impact parameter measured along 
the major axis {p) and the distance along the minor axis 
where the line of sight intersects the disk plane (yo). 


P 


1 




Vo 


\/i + {yo/p)^ 


( 3 ) 


These values can be determined from the impact param¬ 
eter (p) and the azimuthal angle in the plane of the sky 
(a). 

psina 

p = pcosa ; yo = -- (4) 


cosz 


We note that substituting the values from equations [3| 
and 1^ into equation [^yields the expression in equation 
Based on the above, we find Vrotivios = 2.1, p = 
22 kpc, yo = 37 kpc, r = 43 kpc, and (j) = 59°. We 
mark this location with a red circle in Figure which 
shows the projection of the disk on the plane ofthe sky 
along with y—z (edge-on) and x—y (face-on) views of the 
disk plane. The velocity ratio implies deprojected rota¬ 
tion velocities Vrot = [—410, —540, —760] km s“^ for the 
three absorption components. All of these values would 
be inconsistent with the Tully-Fisher r elation (e.g., Tully 


& Fisher 1977 McGaugh et al. 2000) and t he o bserved 
rotational velocity ot the stellar disk (Section 3.2), which 
suggests that the absorbers are either not located in the 
disk plane defined by the stellar disk, or their kinematics 
are not dominated by rotation. 


4.3. Kinematics of gas in a warped disk 

Our line of sight would also intersect clouds located 
above and below the disk plane, and if the kinematics 
were dominated by rotation, the Vrot/vios ratio would 
also be determined by equation]^ The relevant (j) angle 


would depend on the y coordinate along the minor axis, 
which can be calculated for any point along the line-of- 
sight vector L by considering its projection along the y 
and z planes. 


X = p ; y = yo — [LI sinz ; z = zq + \L\ cosz (5) 

Here the x coordinate is measured along the major axis 
in the disk plane, the y coordinate is measured along the 
minor axis in the disk plane, and the z coordinate is mea¬ 
sured along the vertical axis above the disk plane. For 
example, the reference position z/g = 37 kpc corresponds 
to zg = 0 (i.e., crossing the disk plane), and a location 
jL] = 40 kpc along the line of sight from this reference 
position would have y = 0 and z = 16 kpc (i.e., crossing 
the major axis of the disk). The y coordinate then de¬ 
termines Vrot/vios for every point along the line of sight 
based on equations and |3] 

Therefore, we can determine the location that would 
yield an observed line-of-sight velocity for a given rota¬ 
tion velocity. If we use the rotation velocity of the stellar 
disk (vrot = 269 km s“^), the strongest absorber could be 
at jL] = 21 kpc (z = 8 kpc, r = 30 kpc, Vr^vios = 1.4), 
which is marked by a red star in Figure [m If this lo¬ 
cation were associated with a warped disk, the relevant 
warp angle would be tan“^ z/r = 16°. 


4.4. Kinematics associated with outflowing gas 

If the kinematics of the absorbers were dominated by 
an outflowing wind, we could convert the observed line- 
of-sight velocities into three-dimensional wind velocities. 
The ratio of the wind velocity (vw) to the line-of-sight 
velocity (vios) can be expressed in terms of the angle 
(iw) between the wind vector and the line of sight. 


^los COS tuj 

For a wind vector that is orthogonal to the disk plane 
(i.e., along the z axis) the angle iw would be equal to 
the disk inclination angle (z). More generally, the angle 
Zu, can be expressed in terms of the length of the wind 
vector {vru) and its projection on the plane of the sky 
(i.e., the impact parameter, p). 


sinzu, = 


( 7 ) 


Following [Gauthier & Chen| (|2012D , can be expressed 
in terms ot the z height above the disk plane and the 
half opening angle of the wind with respect to the z axis 
(0»). 

rw = ( 8 ) 

COS 0^ 


Thus for any location along the line of sight (equation]^, 
we can calculate the galactocentric radius (r^,), the open¬ 
ing angle of a wind vector that intersects the line of sight 
(0U,), and the associated Vw/vios ratio. 

We can therefore determine the minimum opening an¬ 
gle required to intersect the line of sight and the max¬ 
imum opening angle beyond which one would expect a 
broader range of line-of-sight velocities than is observed. 
The minimum opening angle = 51° corresponds to 
jL] = 55 kpc (z = 21 kpc, Vw/vios = 1-6), and we 
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mark this location with a red triangle in Figure For 
an outflow with an opening angle du, > 51°, the hne of 
sight would enter (exit) the outflow cone at some lower 
(higher) z height for which a smaller (larger) f raction of 
its velocity wo uld be along the line of sight (e.g., [Gauthier 
&: Chen|20^). For example, an outflow with = 52"" 
would span a range |L| = 46-68 kpc {z = 18-26 kpc) 
and exhibit v^/vios values that a vary by a factor of 
1.9, which matches to the range of absorption-line veloc¬ 
ities we observe. Therefore, in the context of a conical 
outflow model, larger opening angles 6^ > 52° can be 
ruled out because they would be associated with an even 
broader range of line-of-sight velocities than we observe. 
We note that an outflow model with > 60° would also 
be associated with redshifted components, and a model 
with 9uj > 67° would additionally produce a blueshifted 
“down-the-barrel” component in self-absorption towards 
the foreground galaxy, neither of which is observed. 

5. DISCUSSION 

We have presented results on a Milky-Way-like disk 
galaxy at z = 0.4, including the column density and kine¬ 
matics of its circumgalactic gas as traced by Mg ii and 
Fe II absorption lines at impact parameter p = 27 kpc. 
The CGM absorption lines are blueshifted in the sense 
of the galaxy’s rotation, and here we evaluate the hy¬ 
pothesis that this gas is associated with an extended, 
warped gaseous disk. We also consider whether some 
components of this absorption could be associated ex- 
traplanar halo gas. This interpretation is relevant be¬ 
cause an extended disk may trace rec ent accretion of gas 
from the intergalactic me dium (e.g., Ostriker & Binney 
1989 Stewart et al.|2013 l, and extraplanar halo gas may 
trace outflows or gas recycling in a galactic fountain (e.g ., 
Shapiro fc Field|fl976 Breginan]|1980| Ford et ar||20I4 ). 

■amts t 


We also discuss the size constraints that can be placed 

on the absorbing clouds due to the extended size of the 
background source. We present this discussion in the 
context of previous observations of gas around galaxies 
as traced by emission and absorption, and we compare to 
expectations from simulations of gas in the circumgalac¬ 
tic medium. 

5.1. Evidence for an extended, warped disk 

In the local universe, extended H i disks have been 
observed in emission around disk galaxies with spatial 
extents rn i = 10-50 kpc for mass surface densities 


Eh I = 1 Mq pc“^ {Nh « 10^° cm“^); this spatial ex¬ 
tent is a function of H i mass such that galaxies with 
Mh t > 10^° Mp) ge nerally have rn i > 30 kpc (e.g., 
Broeils & Rhee 19971. Most of these H i disks exhibit 


In the context of this paper, an extended disk provides 
a straightforward explanation for the optically thick ab¬ 
sorption component at vios = —195 km s“^. Given the 
incidence of extended disks around massive, star-forming 
galaxies in the local universe, it is reasonable to expect 
that our line of sight would intersect an extended disk. 
Galaxies at z = 0.4 also tend t o have higher gas fractions 
than galaxi es at z = 0 (e.g., Rao et al. 2006 Combes 
et al. 20131, which provides additional support tor this 


interpretation. Furthermore, the estimated H i column 
density is consistent with the expected surface density of 
an extended disk, and the kinematics agree with a flat 
rotation curve (vrot ~ 270 km s“^. Section 4.31. The 
velocity dispersion a = hj^fi = 17 km s“^ is consistent 
with the vertical kinematics of a relatively thin disk. 

5.2. Constraints on an outflowing wind 
Galactic winds are known to be common among star¬ 


forming galaxies at z ^ 0.5 (Martin et al. 2012 


Rubin 


eraI]|20T4l ), so it is worthwhile to consider whether the 


absorption we observe could be associated with an out¬ 
flowing wind. Estimates of the half-opening angle for 
winds range fr om 6>„, » 30° for nearby galaxies like M82 
and N GC 253 (Heckman et aHl990l Westmoquette et al. 
,20111 to 0^, « 40°-60° for larger'samples of star-forming 
galaxies at z = 0.1-1 (jChen et al)]|2010c Bordoloi et al. 


2011||K acprzak et al.||2012t [Martin et al. 2012 [Rubin et 

al.||2014p. In addition, Martin et al. (|2012 ) find that the 

opening angle may depend on velocity, such that higher 
velocity outflows tend to be more collimated, while Ru-| 
bin et al. (2014) find evidence that galaxies with larger 
star-formation rate surface densities (Esfr) tend to have 
larger wind opening angl es. 

As discussed in Section |4.4[ only a small range of wind 
opening angles 51° < < 52° could intersect our light 

of sight without producing a broader range of line-of- 
sight velocities than is observed. Furthermore, this would 
require a very fortunate geometry (i.e., the background 
source would have to be almost exactly at the edge of the 
outflow cone), and it would also require a large opening 
angle for a > 300 km s“^ outflow. In particular, the 
minimum opening angle 0^, = 51° would require an out¬ 
flow velocity ~ —320 km s“^(uu, « —580 km s“^) to 
explain the lowest (highest) velocity component we ob¬ 
serve. Furthermore, the maximum opening angle 0^, = 
52° would require Vu, « —460 km s“^ to reproduce both 
the low-velocity and high-velocity components. While 
there is evidence that galaxies with high Esfr valu es can 
produce large-scale, high-velocity outflows (e. 


man et al.||2011[ IDiarnond-Stanic et al.||20i2[ [Law et 


op 


1976 Bosma 1991| jSandsi et al. 2008), and the angles 


Heck- 


^ 2012| [Sell et al.||2014p , the foreground galaxy has 


warp s beyond the optical extent of the disk (e.g., Sancisi Esfr = SFR/fhrri) = 0.01 Mq yr This is an order 


of these warps witn r espect to the stellar disk can be as 
large as 20°-30° (e.g., Garcia-Ruiz et al.|2002 ). In terms 
of kinematics, these ■warped disks are also dynamically 
thin, with typical velocity d ispersions « 10 km s“^(e.g. 


winds (Esfr > 0.1 Mq yr 

Heckman 

2002) 

approximately equal to the sma 

lest Esfr 

value 


z ~ 0.5 galaxies that have been tar geted in absorption- 


Bosma|1991 Bottema|1996 ). There are several explana- 
tions tor the physical origin of warps, including interac- 
tions with nearby dwarf galaxies (e.g., Shang et al.|1998 ) 
and cosmi c infall of gas with misal i gned angular momen¬ 
tum (e.g., Ostrik erfcI^ ney|1989[ Jiang fc Binney|1999 


line studies of gala ctic winds (e.g., Kornei et al. 2012 


Shen fc Sellwood||2006 ). 


Rubin et al.j 2014|). While it is clear that large-scale 
galactic winds are capable of producing strong Mg ii ab¬ 
sorption over a wide range of azimuthal angles around 
star-forming galaxies with high Es fr values (e.g., Lund- 
gren et al. |2012l|Rubin et al.|2014| ^, it is unlikely that the 
absorption lines we observe around this Iow-Esfr galaxy 
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Fig. 5.— Top Left: Visualization of the disk plane projected on the plane of the sky. The solid and dotted lines show a circular r = 43 kpc 
disk in the same plane as the stellar disk. The x and y arrows indicate a coordinate system with x values measured along the major axis 
and y values measured along the minor axis. Top Right: Visualization of a wind model with opening angle 9^ = 51° projected on the 
plane of the sky. The solid and dotted lines correspond to locations in the outflow cone with 2 : heights of above the disk plane of 5 kpc, 
10 kpc, 15 kpc, and 21 kpc. In this model the line of sight would intersect the edge of the outflow cone at 2 = 21 kpc. Bottom Left: A 
face-on view of the disk plane with an r = 10 kpc black circle that represents the stellar disk, an r = 43 kpc black circle that represents an 
extended disk, and a dashed line that represents the line of sight. The red points correspond to three potential scenarios that could explain 
the strongest abs orption component at u = —195 km s“^. The red circle corresponds to an extended disk in the sa me p lane as the stellar 
disk (Section |4.2[ l. The red star corresponds to a warped disk that is consistent with a flat r otation curve (Section |4.3[ l. The red triangle 
corresponds to tne edge of an outflow cone that could intersect the line of sight (Section |4.4[ |. Bottom Right: An edge-on view of the disk 
plane with the same components as the left panel. The r = 10 kpc is drawn somewhat thicker than the r = 43 kpc disk for visualization 
purposes. The 2 height is defined relative to the plane defined by the stellar disk. 


are associated with an outflowing wind. 

5.3. Galactic fountains and recycled accretion 

The fact that we are tracing circumgalactic gas via 
Mg II and Fe ii absorption lines indicates that this gas 
must have been enriched by previous generations of su¬ 
pernovae (i.e., it is not pristine, metal-free gas). It is 
therefore worth considering models of a galactic foun¬ 
tain, which describes the kinematics of supernova- driven 
gas flows that eventually return to the disk (e.g., Breg- 


man|1980 Collins et al.|2002 Fraternali fc Binney|2006 ), 
and models of recycled accretion in cosmological simula- 
tions, which suggest that gas infall at 2 : < 1 is dominated 
by material that was previously ejected in a wind (e.g., 


Oppenheimer et al. 2010 iFord et S(]|2014 ). 


The Bregman 1 1980) galactic fountain model describes 
how hot gas moving vertically upward and radially out¬ 
ward would condense via thermal instabilities into clouds 
that subsequently fall back towards the disk. If no ex¬ 
ternal torques were applied to the gas, the orbital angu- 
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Fig. 6.— Similar to Figure]^ the left and right panels show face-on and edge-on views of the disk plane. The colored points correspond to 
physical scenarios that could explain each of the observed velocity components . Th e red star corresponds to a warped disk, which provides 
the best explanation for the strongest component at u = —195 km s“^ (Section |5.l| . Th e gr een triangle corresponds to co-rotating gas near 
the disk-halo interface that could explain the component at u = —256 km s"'^ (Section |5.5| . The blue square corresponds to the potential 
location of a tidal stream that could explain the component at v = —363 km s"'^ (Section |5.5[ |. The inset in the right panel illustrates 
these three absorption components in velocity space. 


lar momentum would be conserved such that gas clouds 
would follow a ballistic orbit before returning to the disk 
near the radius where the gas originated. In the con¬ 
text of this paper, such a model would not be able to 
explain the large specific angular momentum we observe 
at r « 30 kpc. In other words, the rotational velocity of 
supernova-heated gas originating from r < 5 kpc would 
decrease by a factor of six as it traveled radially outwards 
(e.g., decreasing from Vrot = 270 km s“^ at r = 5 kpc to 
Vrot = 45 km s“^ at r = 30 kpc). In contrast, the rota¬ 
tional velocities we infer at r « 30 kpc are comparable 
to the rotational velocities at r « 5 kpc. Thus it is clear 
that a galactic fountain model in which the angular mo¬ 
mentum of ejected gas is conserved would be inconsistent 
with our observations. 

While the expectations for gas kinematics b ased on 


purely ball i stic orbits are well unders tood (e.g., |Collins 
et al.|[20M Fraternali &: Binney||200(i ), the angul ar mo- 


mentum of gas in a galactic fountain can also be modified 
by hydrodynamic interactions with the gaseous halo. For 
example, several authors have suggested that the drag 
associated with motion through the halo would absorb 
angular momentum from the fountain gas (e.g 


Frater¬ 


nali &: Binney [2008 Melioli et al.||2008 Marinacci et al. 
201ip . This effect would increase the tension between 


model predictions and our observations of high specific 
angular momentum. On the other hand, some models 
have suggested that rapidly spinning halos could im¬ 


part ang ular momentum to ejected gas (e.g.. Brook et 
al. 2012). Furthermore, models of outflows and inflows 


m some cosmological simulations suggest that most of 
the Mg II gas in g alaxy halos is ass ociated with recy¬ 
cled accretion (e.g.. Ford et al.||2014|. In that sense, we 
cannot rule out the possibility that the gas we observe 
at r « 30 kpc could have previously been ejected in a 
wind. That said, it is unclear whether the specific an¬ 
gular momentum of ejected gas could be increased by a 


factor of six to match our observations. A more definitive 
answer would require better understanding of the kine¬ 
matics of hot halo gas and its hydrodynami c interactions 
with ejected gas (e.g., Barnabe et al.|[2006). 


5.4. Comparison to models of cold-mode accretion 
While there is still debate about the fraction of gas 


temperature of the halo (e.g., 

White & Frenk|1991 iKeres 

et al. 

2005 

Nelson et al. 

20 

13), cosmological hydrody- 


namic simulations generically predict that there should 
significant infall of lO'^-lO® K gas onto galaxies, partic¬ 
ularly for halo masses Mh < 10^^ Ad©. The filamentary 
structures that provide gas to Milky Way-mass galaxies 
at z < 1 are expected to be less dense than filaments 
at z > 2 and less likely to survive the journey inwards 
from the virial radius to the outskirts of the galactic disk 
(e.g., Keres et ar||2009 ). Nevertheless, some simulations 
suggest that a traction of the infalling cool gas at z < 1 
may reach the outer disk and bring with it a significant 
amount of angular momentum (e.g. , Rosk ar et al. 2010 
2013 Danovich et al.|[2015). As discusser 


Stewart et al. 


in Section |5.1[ the existence of warped disks in the lo¬ 
cal universe has been interpreted as evidence for cosmic 
infall, so it is worth considering how the kinematics we 
observe compare with predictions for infalling gas from 

zoom-in simulations. _ 

One r elevant comparison is with [Stewart et al. 
( 2011a|b ), who suggest that the end of cold-mode ac- 
cretion at z < 1 for Mh ~ 10^^ Ad© halos could be ob¬ 
served via Mg II absorption systems that co-rotate with 
the galactic disk. When tracking the angular momen¬ 
tum of infalling gas in their simulations, they find that 
recently accreted cool gas has a spin parameter « 4x 
larger than that of the dark matter. In their simula¬ 
tions, the most recently accreted gas has the highest spe¬ 
cific angular momentum and forms an extended, warped 
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disk with velocities that are offset from system by ±100- 
300 km s“^. The gas is expected to sink towards the 
center of t he halo within 1-2 d ynamical times (1-2 Gyr 
1 ; 


Stewart et al. 


2013). While the interactions 


at z 

between outflows and inflows of hot and cool gas are not 
yet well understood, the properties of the cold flow disks 
in these simulations are broadly consistent with our ob¬ 
servations. 

5.5. Comparison with high-velocity clouds, tidal 
streams, and extraplanar gas 

In addition to the dominant component of neutral gas 
in an extended disk, deep H i observations have revealed 
that galaxy halos in the loca l universe often contain 10^- 
10® Mcr, of neutral gas (e.g., Boomsma et al. 12005 Oost- 

i ^ T n Frf • - ri ^ i "-' i - 


erloo et al. 2007 Sancisi et al. 2008 'i'his extrapla¬ 


nar halo gas has been characterized in the Milky Way 
through the study of high-velocity clouds (HVC s; e.g. 


Wakker fc van Woerden|1997 Putman et al. 20121, which 
trace Mh i « 3 x 10® Mq of neutral gas associated with 
the Magellanic Stream, Leading Arm, and Magellanic 
Bridge, plus an additional Mh i » 3 x 10^ Mct) d is- 


The tvoical mass of an HVG com 

plex 

10® Ma> { 

Wakker & van Woerden 

1997 

2012 

1, anc 

cloud complexes as smal 

1 as 3 


Putman et al. 
H I « lO'^ have 


been found in deep H i observation s of M31 (e.g., Thilker 
et al.||2004 Westineier et al.|[2007 |. 


The most prominent gaseous feature in th e halo of the 


10^ Aip) of neutral gas (e.g., Putman et al.|2003 Nidever 
et al.|[Ml0|. The velocity offset between the Magellanic 


System and the M ilky Way (^lmc ~ 380 km s Kalli- 
vayalil et al.|2006 ) is comparable to what we observe tor 
the optically thin component at vios = —363 km s“^. In 
addition to this systemic velocity offset, the Fe ii column 
density and velocity width we find are comparable to the 
observed pro perties of Fe ii absorbers in the Magellanic 
Stream (e.g., Fox et al.|2010 |. In that sense, the data are 
consistent with this high-velocity component being asso¬ 
ciated with a Magellanic Stream-like structure in the 
halo of the foreground galaxy. We visualize the poten¬ 
tial location of such a stream in Figure with a blue 
square that is d = 50 kpc away from the center of the 
foreground galaxy, approximately equal to the distance 
betw een the Milky Way and the Large Magellanic Gloud 
(e.g., Freedman fc Madore||2010 ). 

The extraplanar halo gas observed within Az ~ 10 kpc 
of nearby disk galaxies often shows evidence for co¬ 
rotation with the thin disk, and there is evidence that 
gas at larger z h eights above the disk tends to rotate 
more slowly (e.g. , lRand||1997| |Swaters et al.||1997t |Fra- 
ternali et al.|2002 1. 't his extraplanar component that has 
a connection to disk rotation is som etimes called “disk- 
halo” gas (e.g., Putman et al. 2012). In the context of 
this paper, it is conceivable that the optically thin ab¬ 
sorber at vios = —256 km s“^ could be tracing an extra¬ 
planar gas cloud with kinematics dominated by rotation. 
For example, locations near |L| = 40 kpc (z = 16 kpc, 
r = 27 kpc) have Vrot/vios ± IT, implying that this 
velocity component could have a rotational velocity as 


small as Vrot ~ 270 km s“^, which would be consistent 
with the disk rotational velocity. This location is marked 
with a green triangle in Figure]^ Such a configuration 
would imply essentially no change in Vrot with z height 
up to Az « 8 kpc above the warped disk, which would 
be in contrast with the «15 km s“^ kpc“^ decrease in 
rotation veloci ty with z height for nearby galaxies like 
NGG 891 (e.g., Heald et al.|2006 ). That said, the precise 
location of this cloud along the line of sight cannot be 
uniquely determined from the existing data, so it is not 
possible to draw definitive conclusions about the verti¬ 
cal structure of co-rotating halo gas. Nevertheless, this 
intermediate-velocity absorber does exhibit the general 
characteristics of a co-rotating extraplanar cloud near 
the disk-halo interface. 

5.6. Constraints on cloud sizes 

The fact that we are observing gas in absorption to¬ 
wards an extended background source provides an im¬ 
portant constraint on the transverse size of the absorb¬ 
ing clouds. Previous studies have estimate d clouds sizes 
based on photoionization modeling (e.g., Rigby et al. 


2002 Simcoe et al. 12006 ) and based on differences be¬ 


tw een spectra of multiply lens ed quasars (e.g., Rauch et 


al.||1999| |Petitjean et al. 2000|. These studies often find 
Houd sizes 1 = l\nInn 100 pc (e.g., Schaye et al.|20'07 


thoug 

h sizes estimates as 

et al. 

1999 

Rigby et al. 

10 kpc (e.g., 

Stocke et ai. 


Wrixon 

1972 

Mathewson et al. 

1974 

1, which extends 


more than lUC 

"" across the sky and contains more than 


10 pc (e.g., Rauch 
and as large as t = 1- 


eF^[20l4 | are not uncommon. 


Lehner et al.|2013 Werk 


All three velocity components we observe are well fit 
by models with unity covering fraction, implying that 
the absorbing clouds cover the entire background source. 
This is clear for the strong, saturated Mg ii and Fe ii 
lines at V = —195 km s“^ for which there is no resid¬ 
ual flux at line center. For the weaker components at 
V = —256 km s“^ and v = —363 km s“^, the depth of 
the Mg II A2796 lines (r « 2) approach zero, and the 
relative depth of the Mg ii A2803 lines (t « 1) confirms 
that the Mg ii A2796 lines are not highly saturated. In 
other words, the small residual flux at line center for 
Mg II A2796 in these weaker components is explained by 


a moderate column density (Am 


g II 


10 


13 


cm 


-2 


) with 


unity covering factor rather than a higher column density 
with partial covering of the background source. 

Regarding the size of the background source, we use the 
radius that contains 90% of the total flux (rgn). Based 
on the best-fit de Vaucouleu rs profile (Diamond-Stanic 
et al.|2012||Geach et a l.|20T4t, this corresp onds to rgn = 
520 pc [rgo/r^Q = 5.5; Graham et al.||2005 ). Accounting 
for the difference in angular diameter distance between 
Zbg = 0.712 and Zfg = 0.413, the relevant size scale at 
the redshift of the foreground galaxy is I = 390 pc. Using 
this as an estimate of the size of the absorbing clouds, we 
can calculate the volume density of the gas {uh = Nh/1) 
and the cloud mass (Me = 4/37r(Z/2)^nHmp/i.) from the 
estimated column density. 


nm -sf ^^390 pC^ 

Uh = 0.01 cm -——^ -;— 

V1.2 X 101® cm-2 J \ I J 


Mr — 8 Mlo 


Nh 


® I 1.2 X 101® cm-2 J I 390 


/ I 


pc 


(9) 


( 10 ) 
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The fiducial value Nh = 1-2 x 10^® cm“® is for the 
strongest component, assuming no dust depletion and 
solar metallicity. In comparison, the Nh, nn, and Me 
estimates for the weaker components ar e sm aller by a 
factor of 30-40. As discussed in Section [41} these Nh 
values would be larger if the gas-phase abundance of Fe ii 
were reduced by dust depletion or sub-solar metallicity. 

The most robust constraint on cloud properties asso¬ 
ciated with our analysis is the I = 390 pc lower limit 
on the transverse cloud size. We have assumed in equa¬ 
tions M and m that this value for the transverse size is 
comparable to the absorption path length, which may or 
may not be accurate. Our size estimate is consistent 
with the range of cloud sizes found in previous stud¬ 
ies, although there is some tension with smaller size 
estimates fo und previously fo r weak Mg ii abs orbers 


{1 ~ 10 pc; 

Rauch et al. 

1999 

Rigby et al. 

2002 

). Our 

results suggest that iow-ionizat 

ion clouds with densities 


as small as nn = 3 x 10“^ cm“^ and masses as small as 
Me = 0.2 M.Q can survive in galaxy halos at distances 
r = 30-50 kpc. 


5.7. Implications for the geometry and kinematics of 
circumgalactic gas 

A number of previous studies have combined high- 
resolution spectroscopy of background sources with low- 
resolution spectroscopy of foreground galaxies to test 



y., Steidel et aL|2002 

Kaeprzak et al.|2010l 

Bouche et 

al. 

2013 Burchett et al. 

20B'Chen et al.||20l4D. There 


del et al. 2002 Kaeprzak et al. 
disk models can rarely explain a. 


nents can be as s ociated with an extended disk (e.g., Stei- 

[2M0) 

.1 of thi 


2011b), but pure 
e absorption com¬ 


ponents. This is consistent with earlier work on the kine¬ 
matics of Mg II absorbers, which found that the data 


and weaker “halo” 

components (e.g.. 

Briggs et al. 

1985 

Lanzetta & Bowen| 

1992| |Charlton & Churchill||19ti 

8). 


panding the modeling to include warped disks, which 
are known to be ubiquitous among disk galaxies in the 
local universe. The existence of warps has implications 
for estimates of the galactocentric distance and three- 
dimensional kinematics associated with extended disks, 
and for the range of azimuthal angles over which disk- 
related absorption could be detected. In addition, we 
show how individual absorption components that are not 
associated with an extended disk can be associated with 
other structures that are known to exist in galaxy halos, 
including extraplanar gas and tidal streams. Further¬ 
more, we demonstrate how this science can be done us¬ 
ing background galaxies rather than background quasars, 
which provides complementary information about the 
sizes, densities, and masses of the absorbing components. 


6. SUMMARY 

We have used high-resolution spectroscopy of an ex¬ 
tended background source to study the circumgalactic 
gas associated with a Milky Way-like galaxy at z = 
0.413. We find evidence for an r « 30 kpc warped disk 
that follows the rotation curve of the foreground galaxy, 
consistent with observations of exte nded H i dis ks in the 
local universe (see Figure|4l Sections |4.3l Section|5.1|). In 
the context of theoretical models, this gas could be asso- 
ciated with recycled accretion that has acquired a signif¬ 
icant amount of angu lar m omentum since being ejected 
> 1 Gyr ago (Section and it is also consistent with 
predictions for infallin g coo l gas that co-rotates with the 
galactic disk (Section ]5.4[ ). In principle, a large-scale 
galactic wind with halt-opening angle 0^, « 50° and 
Vui ~ 500 km s“^ could reproduce similar line-of-sight 
kinematics, but the low Esfr for the galaxy and the 
narrow line-of-sight velocity range for the abso rbers sug¬ 


gest that this is unlikely (see Sections 4.4 


5.2). We also 


detect two weak Mg ii absorbers whose kinematics are 
consistent with being associated with a ti dal stream and 
co-rotating extraplanar gas (see Section |5.5[ Figure]^. 
The extended nature of the background source allows us 
to place a lower limit on the transverse size of the ab¬ 
sorbing clouds {I > 0.4 kpc), which is broadly consistent 
previous estimates in the literature. 
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TABLE 1 

Parameters of model fits to circumgalactic absorption lines 
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Note. — Col. (1): Redshift. Col. (2): Velocity [km s ^]. Col. (3): Ion. Col. (4): Log column density [cm ^]. Col. (5): Doppler 
parameter [km s“^]. Col. (6—11): Transitions considered, along with optical depth at line center and rest-frame equivalent width [A] 
for each transition. 
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